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ABSTRACT 

We perform hafo occupatfon distribution (HOD) modeling of the projected two-point correlation 
function (2PCF) of high-red shift (z ~ 1.2) X-ray-bright active galactic nuclei (AGN) in the XMM- 
COSMOS field measured by lAUevato et alj The HOD parameterization is based on low-luminosity 
AGN in cosmological simulations. At the median redshift of z ^ 1.2, we derive a median mass of 
1.02j;o 23 X 10^^ Mq for halos hosting central AGN and an upper limit of 10% on the AGN 
satellite fraction. Our modeling results indicate (at the 2.5a level) that X-ray AGN reside in more 
massive halos compared to more bolometrically luminous, optically-selected quasars at similar redshift. 
The modeling also yields constraints on the duty cycle of the X-ray AGN, and we find that at z 1.2 
the average duration of the X-ray AGN phase is two orders of magnitude longer than that of the 
quasar phase. Our inferred mean occupation function of X-ray AGN is similar to recent empirical 
measurements with a group catalog and suggests that AGN halo occupancy increases with increasing 
halo mass. We project the XMM-COSMOS 2PCF measurements to forecast the required survey 
parameters needed in future AGN clustering studies to enable higher precision HOD constraints and 
determinations of key physical parameters like the satellite fraction and duty cycle. We find that 
N'^/A ^ 5 X 10® deg~^ (with N the number of AGN in a survey area of A deg^) is sufficient to 
constrain the HOD parameters at the 10% level, which is easily achievable by upcoming and proposed 
X-ray surveys. 

Subject headings: dark matter, galaxies: nuclei, large-scale structure of the universe, AGN: general 



1. INTRODUCTION 

Galaxies are believed to reside in dark matter ha- 
los and as such make excellent probes of str i icture 
formation in the universe (e.g.. I White fc FrenkI [l99ll: 



Kauffmann et al. l993'; Navarro et al."1995VMo & Whiti 
1996; Kauffmann et al. 1999; Springel et al. 2005). Ob- 
servations suggest that there exists a supermassive 
black hole at the center of every massive galaxy (e.g., 
iSoltanI I1982D . and hence a physical connection be- 
tween black holes and their host dark matter halos 
is naturally expected. The black hole-halo relation- 
ship has been widely studied in both semi-analytic 
models and cosm ological hydrodyna mic simulations 
(e.g. iKauffmann fc Haehnclt 200 01: iWvithe fc Loebl 
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Conrov fc Whitd |2012[ ). The spatial clustering of 



black holes can be used to probe this relationship and 
to obtain constraints on the formation and evolution 
of supermassive black holes and their role in galaxy 
formation. 

Clustering analyses of different types of active galac- 
tic nuclei (A GN) have been wid e ly studied in the lit- 
erature (e.g., iCroom et all 120041 : iPorciani et all 120041: 

I j onat han. richardson@uchicago . edu | 



I Croom et all 120051 IGilli et al.l 120051: iMvers et al.l mim 



2007a, b; Coil et al. 2007; Shen et al. 2007; Wake et al 
2008; Shen et al. 2009: Ross et al. 2009: Coil et al 
20091: iHickox'et all 2009, 2011; AUevato et al. 12011 



Donoso et al.ll201^K rumpe et al..i201Q:,Cappelluti et alj 
2012; Whi te et al.l l2012: SheneLalJ|20]J) . Many of these 
studies employ samples of optically-bright quasars, which 
provide large sample sizes and accurate redshifts that 
enable precise measurements of the relationship between 
the growth of black holes and their parent dark mat- 
ter halos over a wide range in redshift. However, while 
quasar clustering has provided a wealth of information 
on the cosmological evolution of black holes, samples of 
optical quasars preferentially probe massive black holes 
accreting at high Eddington rates (Hopkins et al. 2 0091 ; 
iKellv et all 120 101 ). To fully understand the growth of 
black holes over cosmic history, we require a systematic 
study of the relationship between black holes and their 
host halos spanning a larger range of redshifts, luminosi- 
ties, and black hole masses. This can be achieved through 
X-ray surveys. 

AGN clustering is most often studie d using th e two- 
point correlation function (2PCF; e.g.. IArDlll970f ). For 
a given cosmological model, the amplitude of the large- 
scale AGN 2PCF relative to the dark matter 2PCF (i.e., 
the square of the AGN bias) can be estimated. From the 
value of this bias factor an estimate of the typical mass 
of AGN-hosting dark matter halos can then be obtained 
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(;e.g.. lJineill998l:lSheth et~alll2001[ ). Previous studies sug- 
gest that, compared to optical quasars, X-ray-selected 
AGN are more strong ly clustered and reside in more ma- 
sive host halos 
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, ios (e . g.. lUUii et ai.lizuu5l: 
Hickox et al.]l2009l : lAllevato et al.ll20'Tll : 
201211 . However, the host halo mass range is insufficiently 
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constrain ed for a definitive con clusion (for a recent re- 
view, see iCappelluti eTaI|[20T2t ). 

Beyond simple bias measurements, a powerful ana- 
lytic technique known as the halo occupation distri- 
bution (H OD) formalism ( e.g., Ma & Frv 2000: Sclial 
20M iBerlind fc Weinberel [20021 : fZheng et all l200" 



Zheng fc Weinbereil2007( r can be used to investigate the 
full relationship between AGN and their host halos. The 
HOD is characterized by the probability P{N\M) that 
a halo of mass M contains N AGN of a given type, to- 
gether with the spatial and velocity distributions of AGN 
inside individual halos. As long as the HOD, at a fixed 
halo mass, is statistically i ndependent of th e large scale 
environments of halos (e.g.. lBond et al.|[i99lD . it provides 
a complete description of the relation between AGN and 
halos, allowing the calculation of any clustering statis- 
tic (e.g., N-point correlation function, void probability 
distribution, pairwise velocity distribution) at all scales 
for a given cosmological model. This implies that if the 
HOD can be constrained empirically, it will encode all of 
the information that measured clustering properties can 
convey about AGN formation and evolution. 

The HOD approach has been used by several au- 
thors to interpre t AGN and quasar c lustering mea- 
surements (e.g.. iWake et al.l [2 008: Sh en et alj 120101: 
Mivaii et al.* '2 0111: iStarikova et al, ,20 lit lAllevato et al.l 
2011; Richard soTet al.l 120121: IKavo fc Oguril I2012D. hi 
particular Richar dson et al r (|2012D . and IKavo fc Ogurl 
( |2C)12 ) made clustering measurements of optically lumi- 
nous quasars for both small (< 1 Mpc) and large scales 
and performed HOD modeling of the 2PCF to infer the 
relation between quasars and their host dark matter ha- 
los at z ~ 1.4. Similar analysis was carried out by 
IShen et al.l ([20T2h for the quasar-galaxy cross-correlation 
function (CCF). A natural extension of work on the HOD 
of optically luminous quasars is to study the black hole- 
halo relationship for lower luminosity AGN, selected from 
optical or X-ray survey s. In the curren t work , we use the 
2PCF measurements of lAUevato et al.l ()2011|) and a the- 
oretically motivated m odel of the AGN HOD similar to 
that used by Richards on et al.l (|2012( ) to interpret the 
spatial clustering of X-ray-selected AGN at z ^ 1.2. We 
compare our results with previou s studies of the quasar 
HOD by [Richardson et al.l ([2012D and construct an evo- 
lutionary picture of AGN growth and the role of AGN in 
large -scale structure formation (similarly to Hicko x^t al.l 
[20091 but now from the HOD perspective, and at a higher 
redshift of z ~ 1). 

Our paper is organized as follows. In §2, we briefiy 
describe our data sets, the parameterization of the AGN 
HOD, and the theoretical modeling of the 2PCF. We 
present the results of our HOD modeling in §3. Finally, 
we discuss the implications of our results in §4 and sum- 
marize them in §5. Throughout the paper we assume a 
spatially fiat, ACDM cosmology: fJ„ = 0.26, = 0.74, 
= 0.0435, Us = 0.96, erg = 0.78, and h = 0.71. We 
quote all distances in comoving Mpc and masses in 
units of Alpi unless otherwise stated. 



2. DATASETS AND METHODOLOGY 

In this section we briefly describe our clustering data 
set and methodology for the HOD analysis. 

2.1. Clustering Data Set 

We model t he p rojected 2PCF measurements of 
lAllevato et all ([20Tl[ ) for AGN in the XMM-COSMOS 
survey. COSMOS covers 1.4 x 1.4 deg^ equatorial fielcQ 
and uses multiwavelength data from X-ray to radio 
bands. XMM-Newton surveyed 2.13 deg^ of the COS- 
MOS field in the 0.5-10 kev energy band for a total of 
1.55 Ms. The resulting XMM point source catalog con- 
tains 1822 objects, 1465 of which are XMM-COSMOS 
AGN detected in the soft X-ray band (jHasinger et al.l 
[20071 : [Caroelluti et al.l [20071 [2009^ . 

Using the 780 XMM-COSM OS AGN t hat h ave spec- 
troscopic confirmations, lAllevato et al.l (|201lD applied 
a magnitude cut of Iab < 23 to obtain a clustering 
sample of 593 X-ray-selected AGN. The sample spans 
the redshift range 0.0-4.0 with a median redshift of 
1.2. The X-ray luminosities of the AGN span the 
range ^ lO'^^-lO*'^ erg s~^ with a media n luminosity of 
~_6x lO'*'^ erg s~^. We refer the reader to lAllevato et al.l 
(|201lD for a detailed discussion of the 2PCF calculation. 

Correcting for spectroscopic completeness, we find the 
redshift-averaged observational number density of the 
AGN to be (6.7 ± 0.7) x 10"^ Mpc~^, where the 
^ 10% uncertainty corresponds to the characteristic level 
of scatter we observe between the redshift averages com- 
puted by discrete summation and by fitting a decaying 
expone ntial to the observed redshift distribution (see Fig- 
ure 1 of lAUevato et al.ll2011l for the redshift distribution) . 
The 2.13 deg survey area probes a total volume of 0.09 
Gpc^. For this work, we compute the source density by 
averaging over the volume-weighted number density for 
each redshift bin (i.e., the number of sources within the 
redshift range divided by the volume of the redshift slice) . 

2.2. HOD Modeling 

We ad opt the theoretical HO D model for AGN pro- 
posed bv lChatteriee et all (|2012D . which was inferred for 
lower luminosity AGN (Lboi < 10*^^ erg s~^) in c osmo - 
logical simulations. However, [Richardson et all (|2012[ ) 
have found this model to also provide an excellent in- 
terpretation of the clustering of quasars observed in the 
Sloan Digital Sky Survey (SDSS). We therefore apply this 
particular model to moderate-luminosity X-ray AGN as- 
suming a universality in the general HOD properties of 
AGN. 

This model represents the AGN mean occupation func- 
tion as the sum of its physically illustrative central and 
satellite components, {Nccn{M)) and {Nsat{M)), respec- 
tively. The mean occupation function is given as the 
sum of a softened step function (central component) and 
a rolling-off power law (satellite component). 
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where A/min, ciogAf, Mi, a and Mcut are free parame- 
ters. In order to minimize the parameter degeneracy in 
our modehng, we fix Mcut to a value ^ 10^^ Mq so 
that the exponential factor, exp {— Mcut / M) , approaches 
unity over the range of physical halo masses. This re- 
duces our satellite parameterization to a two-parameter 
power law. 

Using the routine developed in lZheng et al.l (|2007D . we 
perform a Markov Chain Monte Carlo (MCMC) model- 
ing of the projected 2PCF to sample the four-dimensional 
parameter space of the AGN HOD. Each MCMC chain 
discussed herein contains 60,000 points in the HOD pa- 
rameter space. Since the clustering sample is sparse, 
we calculate the value of each point using only the 
diagonal elements of the covariance matrix (see, e.g ., 
the appendixes of IMvers et' al. 2007al: IRoss et al.l 120091 ). 
We fit both the observed 2PCF and number density of 
AGN to constrain the HOD. This routine includes the 
effects of halo exclusion, nonlinear clustering and scale- 
dependent halo bias, where halos are defined as objects 
with a mean density of 200 times that of the background 
universe. The halo mass function is computed accord- 
ing to the formula in i Jenkins et al.l (UOOl) and the large- 
scale halo bi a s fact or is computed using the formula in 
[Tinker et al] (|2005l ). We adopt flat priors in logarith- 
mic space for Mmin and Mi and in linear space for a 
and (TiogM (ciogM > 0). For improved computational 
efficienc y, we further requ ir e tha t 0.5 < a < 4.0. 

As in Richardso n et all (|2012D . we adopt the follow- 
ing conventions in our modeling. We assume that the 
halo occupations of central and satellite AGN are inde- 
pendent (i.e., the number of satellites in a given halo 
does not depend on wh ether there is a central AGN), as 
IChatteriee et al.l (j2012D found no evidence of a correla- 
tion between the activity of central and satellite black 
holes in a hydrodynamic simulation. For halos of a given 
mass, we assume that the central occupation numbers 
obey a nearest integer distribution (|Berlind fc Weinberg! 
120021) and the satellite occupation numbers obey a Pois- 
son distribution. This has also been f ound to be the case 
for AGN in cosmologic al simulations (jDegraf et al.ll20l"TI : 
IChatteriee et al.|[2012l) . We represent the spatial distri- 
bution of satellite AG N within halos as an NFW profile 
(jNavarro et al.lll997D w i th the concentration-mass rela- 
tion from lBullock et all (|200H) . 

where M^ is the nonlinear mass for collapse at z = 0, and 
P = —0.13. We adopt cq = 32, motivated by the high 
conce ntration observed for local AGN (e.g.. ILin fc Mohil 
I2007h . We have verified that our modeling only weakly 
depen ds on Cn for a wide ran ge of values, from ^ 10 to 
~ 60 (jRichardson et al.ll20Tl) . 

To enhance its statistical power, the clustering sample 
has been constructed over a broad redshift and luminos- 
ity range to maximize the volume and the number of 
sources. The 2PCF obtained in this way can be inter- 
preted as an average over the redshift and luminosity 
intervals. However, extending this interpretation to the 
HOD could be problematic since the modeling uses halo 
properties (e.g., mass function, bias factor) at the me- 
dian redshift. Redshift evolution of the halo properties 



can lead to systematic effects that exceed the statistical 
uncertainties refiected in the modeling results. 

Although a full interpretation of the measured 2PCF 
would require additional param eters accounting for evo- 
lution, [Richardson et all (j2012D have demonstrated that 
the HOD can be meaningfully interpreted as represent- 
ing the HOD for objects at the median redshift of the 
sample (to within the quoted uncertainties), if the 2PCF 
measured over a wider redshift range around the median 
redshift (i.e., the average 2PCF) is statistically consistent 
with the true 2PCF of objects at the median redshift. 
Since we are not in a position t o exp licitly check this 
condition for the lAUevato et al.l (j2011l) clustering sam- 
ple, we adopt the median-redshift interpretation with the 
assumption that t he clustering ev olves weakly with red- 
shift. See §4.1 and lRichardson e~ al. (2012) for additional 
discussion regarding this chosen interpretation. 

3. RESULTS 

We first present the constraints on the AGN HOD ob- 
tainable with current data. Then we forecast improve- 
ments in the precision of the 2PCF to estimate the req- 
uisite observational parameters to obtain high-precision 
(few percent-level) HOD constraints. 

3.1. Current Constraints 

In Figure [1] we show our four-parameter HOD fit for 
X-ray AGN at z ~ 1. 2. Panel (a) shows th e projected 
2PCF measurements of lAUevato et al.l (j2011| ) against the 
theoretical prediction of our best-fit HOD model (dot- 
ted line). We identify the best-fit model as the point 
in our four-dimensional parameter space associated with 
the global minimum. If we rank all of the models 
by their statistic in ascending order, the first 68% 
of the models give the range of predicted Wp indicated 
by the shaded envelope. Our HOD model reproduces 
the clustering with a reduced statistic of 0.34 (total 
= 2.39). Given a distribution for seven degrees of 
freedom, the probability of randomly drawing a value 
less than or equal to 2.39 is only 0.065. This indicates 
that while the model successfully reproduces the data, 
the uncertainties quoted on the 2PCF are likely conser- 
vative. We discuss this issue further in H4.ll 

Panel (b) shows the mean occupation function from the 
best-fit HOD model, decomposed into its central (dashed 
line) and satellite (dot-dashed line) components. Simi- 
larly to panel (a), the shaded regions indicate the range of 
the mean occupation function given by the 68% of mod- 
els with the smallest statistic. The panel shows that 
while reasonable constraints can already be obtained on 
the central occupation with current data, only an upper 
limit can presently be obtained on the satellite occupa- 
tion. This upper limit suggests that, at low redshift, typ- 
ically only the most massive halos (> 10^'* Mq) host 
multiple AGN. It also yields an upper limit on the satel- 
lite fraction of X-ray AGN to be < 0.1, where the satellite 
fraction is defined as the ratio of the number density of 
satellite AGN (integrated over all halo masses) to the to- 
tal number density of all AGN. However, we note that 
this interpretation is sensitive to our choice of satellite 
parameterization. Al ternative models of the mean occu- 
pation function (e.g.. lKavo fc Ogurill2012D allow multiple 
AGN to reside in low-mass halos. 
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Fig. 1.— HOD fit to the 2PCF of z ~ 1.2 X-ray AGN. Panel (a); the projected 2PCF of lAUevato et al.l l(20Tll : data points and error 
bars) against the prediction of our best-fit HOD model (dotted line). Panel (b): the mean occupation function of AGN, decomposed into 
its central (dashed line) and satellite (dot-dashed line) components. In both panels (a) and (b), the shaded envelopes indicate the 68% 
confidence regions. Panel (c): the full probability density function of host halo masses for central AGN. This distribution is obtained 
by multiplying the cen tral mean occupation function with the differential halo mass function and averaging over all the models in the 
MCMC chain (see i|3.1l for discussion). Panel (d): the probability density function of the median mass of halos hosting central AGN. This 
distribution is obtained from the median host halo mass of each model in the MCMC chain. The vertical lines denote the central 68% 
confidence interval. 



Panel (c) shows the full host halo mass distribution 
(averaged over all points in the MCMC chain) for cen- 
tral AGN. For each model in the MCMC chain, the dis- 
tribution is derived by multiplying the mean occupation 
function of central AGN, {Nccn{M)) , with the differential 
halo mass function. For a randomly chosen central AGN, 
the curve is the probability density function for the mass 
of its host halo. Finally, panel (d) shows the probabil- 
ity density function for the median mass of halos hosting 
central AGN. At the 68% confidence level, denoted by the 
dotted vertical lines, we find the median host halo mass 
for central AGN to be 1.02tg;^i x lO^^ Mq. We 
do not depict the corresponding distributions for satel- 
lite AGN in this figure — the satellite mean occupation is 
insufhciently constrained for its halo mass distribution 
to be robustly determined. We will discuss these HOD 
constraints further in i)4.2l 

The central mean occupation function shown in panel 
(b) of Figure [T] can also be interpreted as the halo mass- 
dependent duty cycle of AGN (i.e., the fraction of halos 
hosting a central AGN). We estimate the duty cycle fo 
of central AGN around the median host halo mass. For 
each model in the MCMC chain, we average the predicted 



{Ncen{M)) over the mass interval spanned by the central 
68% of the median host halo mass distribution (shown in 
panel d of Figure[T|). From the normalized distribution of 
these mass-averages we infer, at the 68% confidence level, 
an average duty cycle of fo = 0.12 ± 0.02 for z ~ 1.2 X- 
ray AGN. 

3.2. Error Forecasting 

Although the physical constraints obtained in t )3.1l are 
consistent with existing estimates, the current precision 
of 2PCF measurements does not strongly constrain the 
full range of possible HOD models. 

In order to forecast the improvement in HOD mod- 
eling precision that can be achieved with planned X- 
ray surveys, we construct and model a set of simulated 
2 PCF measurements ba sed on the original measurements 
of lAllevato et al.l (|2011D . but with appropriate rescaling 
of the error bars. Since the observational AGN num- 
ber density varies widely between different proposed sur- 
veys, an appropriate Swp scaling relation must account 
for both survey coverage and depth. As will be discussed 
in i )3.3[ for example, a shallow, large-area survey can have 
a surface density (and thus volume density) an order of 
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Fig. 2. — The forecasted HOD modeling precision obtainable as a function of the survey N'^ /A ratio, where is the number of sources 
observed in the survey and A is the survey area. The left panel shows the average fractional uncertainties (see the definition in the text) on 
the 2PCF (dashed line), central mean occupation function (solid line), and satellite mean occupation function (dot-dashed line). The right 
panel shows the fractional uncertainties on physical parameters derived from the mean occupation function; the median host halo mass 
for central (dashed line) and satellite (solid line) AGN and the AGN satellite fraction (dot-dashed line). In both panels, the vertical lines 
denote the N'^ /A value of a particular proposed X-ray AGN survey, as labeled (see Table [T]l. The "XLL" and "CMS" abbreviations refer 
to the XMM-XLL and Chandra 10 deg'^ surveys, respectively. Our erro r fore casti ng i ndicates that a large improvement in HOD modeling 
precision can be realized by planned AGN surveys of modest size. See ^3.21 and i|3.3l for additional discussions regarding the feasibility of 
different planned surveys. 



magnitude smaller than those of deeper smveys. This 
will significantly impact the observed number of AGN 
pairs, which determines the overall precision of the 2PCF 
measurements {5wp oc l/y^Vp^h^, with iVpairs the num- 
ber of pairs). If we assume the redshift distribution to 
be similar for every proposed survey, the number (vol- 
ume) density of a given survey scales approximately as 
the number of sources per unit area. Assuming popula- 
tions with the same clustering, the number of AGN pairs 
would then scale as the product of the total number of 
AGN in the survey, N , and the observed AGN surface 
density, N/A, where A is the angular sky coverage. 

We thus rescale our 2PCF error bars according to 
1/ ^JTP~[A to forecast future HOD constraints. For each 
Tp bin, we take the corresponding Wp value to be the 
value of the best-fit 2PCF shown in panel (a) of Figure [1] 
(i.e., the predicted 2PCF of the HOD model best fitting 
the[Allevato et ah 2011 data). We use the best-fit 2PCF 
values instead of the original measurements to obtain 
a "smoothed" data set free of outliers, since artificially 
reducing the error bar on an outlying data point can po- 
tentially cause the HOD modeling to fail. We fi nd that 
the pa ir separation range probed by the Alievat o et al.l 
(|2011l) sample, 0.5 h'^ Mpc < Tp < 30 h''^ Mpc, is insuf- 
ficient to tightly constrain the satellite occupation, irre- 
spective of the measurement precision. While the central 
occupation is determined by the two-halo term, whose 
contribution dominates the 2PCF above ^ 1 Mpc, 
the satellite occupation is determined by the one-halo 
te rm, which domin a tes at smaller scales. The results 
of iRichardson et al.l H2012D suggest that 2PCF measure- 



ments to at least ^ 0.01 h^^ Mpc scales are required to 
fully constrain the satellite occupation, and we note that 
larger samples will produce more small-scale pairs, al- 
lowing measurement of the 2PCF to smalle r scales than 
was possible with the lAllevato et al.l (|2011[ ) sample. 

Hence, in order to isolate the uncertainty in the satel- 
lite occupation due to measurement precision, we ex- 
tend our dummy 2PCF to smaller scales by adding radial 
bins centered on 0.1, 0.05, and 0.01 Mpc. As with 
the larger radial bins, we assign to each new bin the 
Wr> value pred i cted b y the HOD model best fitting the 
lAllevato et al.l ()2011l ) data. For a small-scale bin cen- 
tered on Tp with upper and lower boundaries at rp2 and 
Tpi , respectively, the uncertainty scales approximately as 

5wp oc Wpl \jT^{rl^ - ?'^i)(27r,„ax + Wp), (3) 

where TTmax is the maximum line-of-sight separation be- 
tween AGN pairs. The constant of proportionality can 
be determined by ma tching this relation to the existing 
lAllevato et al.l 1)20111 ) measurements at Tp < 1 Mpc. 

We now define a statistic to describe the overall preci- 
sion of a given data set in a single number, to which we 
will refer as the "average fractional uncertainty." Ap- 
plied to the projected 2PCF measurements, the aver- 
age fractional uncertainty is calculated by computing 
the fractional error on each data point, 5wp{rp)/wp{rp), 
and then averaging the fractional errors over all rp 
bins. This statistic is calculated similarly for the cen- 
tral and satellite mean occupation functions, {Nccn{M)) 
and {Nsat{M)), respectively, but the halo mass range is 
restricted to the interval over which the mean occupation 
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function is greater than 10~^. We impose this condition 
so as to consider the uncertainty only over the physical 
halo mass range (i.e., the range that contributes signifi- 
cantly to the AGN number density). 

Figure [2] shows the forecasted HOD modeling precision 
obtainable as a function of the survey N"^ /A ratio. The 
left panel shows the average fractional uncertainties on 
the 2PCF (dashed line), central mean occupation func- 
tion (solid line), and satellite mean occupation function 
(dot-dashed line). The right panel shows the fractional 
uncertainties on physical parameters derived from the 
mean occupation function: the median host halo mass 
for central (dashed line) and satellite (solid line) AGN 
and the AGN satellite fraction (dot-dashed line). In 
both panels, the vertical lines denote the N'^ /A value of a 
particular proposed X-ray AGN survey, as labeled. Our 
error forecasting indicates that a large improvement in 
HOD modeling precision can be realized by planned AGN 
surveys of modest size. For a survey of area equal to that 
of the XMM-COSMOS field (2.13 deg^), for example, we 
find that a sample of ~ 1,000 AGN with spectroscopic 
redshifts will be sufficient to double the precision of cur- 
rent constraints, and a sample size larger than 10,000 will 
allow determination of the derived physical parameters 
to to better than 10% statistical uncertainty. In the next 
section we discuss the plausibility of current and future 
X-ray surveys in achieving this precision. 

3.3. Survey Design 

Based on our error forecasts we now evaluate our pre- 
dictions in light of ongoing and future surveys of X-ray 
AGN. In Table [T] we show the relevant source statistics. 
The number counts are obtained usin g the predic tions of 
the cosmic X-ray background model of lGilli et al.l ([2007)0 
and sensitivity curves for various X-ray surveys. All X- 
ray luminosities and fluxes are quoted in the 0.5-2 keV 
band. The counts are for sources detected at all redshifts 
(the redshift distribution peaks at z ~ 1 independent of 
fiux limit). The parameters /20 and /go represent the 
fiux limits at 20% and 80% of the total survey area, re- 
spectively. These two numbers are useful in providing a 
more accurate shape of the sensitivity curve. In column 8 
of Table [T] we compute the approximate source densities 
(number of yield per unit area) of relevant X-ray surveys. 
It is important to note that the larger area, shallower 
surveys have smaller surface densities (and thus volume 
densities) than the deeper surveys. This potentially af- 
fects the number of pair counts in clustering studies. 

We find that the Chandra COSMOS Legacy Survey, 
the XMM-XXL Survey, and a Chandra medium-area 
(~ 10 deg^) survey can all obtain high-precision HOD 
constraints. However, as mentioned previously, our re- 
quirement that spectroscopic redshifts be obtained for 
our sources may prove challenging for a 50 deg^ sur- 
vey like XMM-XXL. Hence we emphasize that a deeper 
survey covering a smaller area is preferable for HOD 
studies, given that the surface density of sources is 
larger in this case and follow-up spectroscopy is bet- 
ter facilitatecQ. It is evident from Table [T] that either 
the COSMOS Legacy Survey or a medium-area Chan- 
dra survey can be utilized to study luminosity evolution 

^ http://www.bo.astro.it/ gilli/counts.html 



of the HOD. The lAllevato et"all (|2011[ ) clustering sam- 
ple has been constructed over a wide luminosity range 
(10*^ — 10^^ erg s~^) and our HOD precision estimates in 
Figure [2] are obtained using the median- redshift /median- 
luminosity interpretation. While these precision fore- 
casts are robust in terms of statistical uncertainties, with 
increasing statistical precision in the 2PCF, systematic 
effects from luminosity and/or redshift evolution will 
become increasingly important and may limit the ul- 
timate precision of constraints obtainable on the HOD 
(i.e., the median- redshift /median-luminosity interpreta- 
tion may no longer apply; see §4.1|) . In this case, more 
sophisticated HOD models with parameters accounting 
for redshift and/or luminosity evolution would be war- 
ranted. 

In both panels of Figure [21 all of the fractional pre- 
cision curves for the HOD modeling scale (roughly) lin- 
early with the precision curve for the 2PCF before fiat- 
tening once it surpasses ~ 10% precision (near N"^ /A ^ 
2 X 10^ deg^). The characteristic shape of these curves 
arises from a transition in the dominant source of statis- 
tical uncertainty in the modeling. As the fractional un- 
certainty of the 2PCF decreases with increasing N"^ /A, 
it eventually falls below the fractional uncertainty of the 
AGN number density, which we have fixed to 10%. We 
investigated a wide range of number densities and frac- 
tional uncertainties (up to a factor of ~ 5 variation) to 
determine their impact on our HOD modeling results. In 
all cases, we recovered similar curves to those shown in 
Figure [2] which flattened near the N'^ /A value of 2PCF- 
number density fractional uncertainty equality. Since the 
planned surveys discussed herein all fall within the num- 
ber density-limited regime {N'^/A > 2 x 10^ deg^), we 
note that their forecasted precisions are limited by the 
assumed 10% uncertainty in the number density, rather 
than the precision of the 2PCF, and therefore may be 
conservative estimates. If the systematics of the surveys 
can be controlled to better than 10% (which we do not 
assume) , it may be possible to obtain slightly tigher con- 
straints on the HOD than those indicated here. 

Finally, we note that the survey error forecasts pre- 
sented herein are calculated for AGN attio-correlation 
measurements. Although it is beyond the scope of the 
current work, some studies have suggested that even 
higher precision clustering constraints may be achievable 
throug h the cross-c orrelation between A GN and galax- 
ies (see lHickox et al, 2009 and lCoil et al.ir2009. for exam- 
ples using X-ray AGN) and could possibly lead to HOD 
constrain ts if the galaxy population is well-ch aracterized 
(see, e.g.. IWake et alll200l IShen et al.llMl ). 

4. DISCUSSION 

We now discuss the systematic issues and theoretical 
aspects of our analysis. 

4.1. Systematic Uncertainties 

As mentioned previously, we apply the median- 
redshift/median- luminosity interpretation of 

Novel statistical techniques have been developed to maximize 
the precision of clustering measurements using photometric redshift 
samples (Myers et al. 2009; Hickox et al. 2011), but the statistical 
power of these measurements is necessarily limited compared to 
spectroscopic studies using comparable sample sizes. 
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TABLE 1 

AGN SOURCE COUNTS IN X-RAY SURVEYS 





Area (A) 


/20 




/so 




Numbcr(Af)/ log[Lx 


/(erg s ^)] bin 


N/A 


N'^/A 




Survey 


(dcg2) 


(erg cm~-^ 


s-i) 


(erg cm~^ 


s-i) 


42-43 


43-44 


44-45 


(deg-i) 


(deg-2) 


Ref 


Chandra COSMOS Legacy 


2 


2.0 X 10- 


-16 


3.7 X 10" 


-16 


1600 


2100 


460 


2000 


8 X 10** 


(1) 


Chandra 10 dcg^ survey 


9 


6.1 X 10- 


-16 


7.9 X IQ- 


-16 


2900 


6300 


1900 


1200 


1.3 X lO'^ 


(2) 


XMM-XXL 


50 


3.1 X 10" 


-15 


1.0 X 10" 


-14 


1300 


5800 


4000 


440 


9.7 X lO** 


(3) 


eROSITA (Deep) 


100 


3.1 X 10- 


-15 


3.2 X 10- 


-15 


4700 


21000 


12000 


380 


1.5 X lO'^ 


(4) 


eROSITA (AH Sky) 


41253 


10" 




10-16 




320000 


1600000 


2000000 


95 


3.7 X 10* 


(5) 



Note. — The values in columns N/ A and N^/A are calc ulat ed over all luminos i ty bins (42 < log[ Lx/(erg s ^)] < 45). Re ferences: 1. F. 
Civan o (private comm. ), 2. Scaled from IKenter et al.l l l2005f) . 3. lElviv et al.l (201% . IPierrel 1 12013) 4. IMerloni et all ll20Ta 'l. 5. IMerloni et all 



IRichardson et alj (|2012l) to our HOD modeling. 
Within this interpretation, the HOD can be taken to 
represent AGN at the median sample redshift if the 
luminosity and redshift evolution of the 2PCF are 
weak (i.e., their systematic effects do not exceed the 
statistical uncerta i nty q uoted on the modeling results). 
IRichardson et all (|2012| ) provide a detailed discussion 
of how intrinsic luminosity or redshift evolution of the 
clustering would affect the HOD modeling. 

To high precision, the clustering of optically-selected 
quasars has not been found to exhibit s i gnificant lu- 
minos i ty evolution (e.g.. Gr oom et al.ll2005t iMvers et ahl 
I2007al: iShen et al.ir2009, .20121 ) . Available studies on the 
luminosity evolution of X-ray AGN clustering do not 
show any evide nce for str ong evolution a t z ^ 1 (e.g., 
iGoil et al.|[2009l ). although iKrumpe et al.l (|2010[ ) demon- 
strated significant luminosity evolution of clustering at 
z ^ 0.25. Future studies with larger samples of X-ray- 
selected AGN are needed to confirm the degree of lumi- 
nosity evolution in X-ray AGN samples. Since the lumi- 
nosity evolution at z '--^ 1 has not been well quantified, 
we assume that the luminosity evolution is weak and is 
incorporated within the uncertainties of the 2PCF. 

There has not been a comprehensive study of the red- 
shift evolution of the clustering of X-ray-selected AGN 
(due to small sample sizes). How ever, we examine the 
redshift evolution of bias from the lAllevato et al.l ()2011[ ) 
sample to infer the approxim ate redshift evolution of the 
2PCF. lAllevato et"all ()20lH) measure the bias of the full 
sample (median (z) = 1.2) to be 6 = 2.98 ± 0.13. They 
then construct a sub-sample spanning a smaller redshift 
interval, but with a similar median redshift {{z) = 1.3) to 
that of the full sample, and evaluate the bias of the sub- 
sample to be 6 = 3.10 ± 0.18. The consistency between 
the clustering amplitudes over the whole redshift range 
and in the narrow redshift range around the median red- 
shift therefore lends supports to our interpretation of the 
H OD modeling resu l ts. 

IRichardson et al.l (|2012D provide a thorough discus- 
sion on the implications of alternative HOD models 
on the quasar 2PCF. In short, IRichardson et al.l ()2012f ) 
found the central occupation to be robust to the chosen 
HOD parameterization but found the satellite occupation 
to exhibit statistically significant model-dependences at 
high halo masses. This systematic effect led to a sta- 
tistically significant difference in some physical param- 
eters. For example, the qua sar satellite fractions ob- 
tained from HOD modeling bv lKavo &: Oguril ()2012D and 



IRichardson et all (|2012[ ) disagreed with each other by 
a large factor (two substantially different HOD models 
were assumed by the two different groups). Although 
breaking this degeneracy in the high-mass satellite oc- 
cupation will require larger samples of small-separation 
pairs coupled with additional constraints (e.g., the distri- 
bution of line-of-sight pair velocities), it is not of central 
concern to the current X-ray AGN constraints presented 
herein. Our modeling is unable to impose strong con- 
straints on the form of the satellite occupation due to 
the lack of small-scale clustering data (< 1 h^^ Mpc). 

As an additional check, we generalize the form of our 
central HOD parameterization by allowing the softened 
step function to approach an arbitrary asymptotic value 
less than unity (see Eg. 1) and repeat the modeling of 
the lAllevato etall ()201lD data. We find that the HOD 
fit is statistically insensitive to this new degree of free- 
dom, and we recover similar relative uncertainties to 
those yielded by the original model. Although our chosen 
model has been calibrated for low-luminosity AGN, we 
emphasize that the current model has successfully repro- 
duced quasar clustering and has been favored by direct 
measu rements of the mean occupation function of X-ray 
AGN (lAllevato et al.llMl) . We compare our HOD con- 
straints obtained from modeling t he 2PCF to those di - 
rectly inferred for X-ray AGN by lAllevato et"aI1 H2012D 
in fjiH 

Finally, while our HOD model successfully reproduces 
the observed clustering of X-ray AGN, it fits the data 
with a reduced statistic of 0.34. This suggests ei- 
ther that our HOD model admits too many free param- 
eters, thereby overly fitting the details in the 2PCF, or 
that the error bars quoted on the 2PGF overestimate the 
true Icr uncertainty. To test whether the four-parameter 
model overfits the data, we fit only the central occupa- 
tion function to the two-halo term (all data points at 
rp > 1 Mpc) with the parameter ciogM fixed to the 
best-fit value yielded by the full modeling. Hence we only 
fit a single free parameter, Mmin, to nine data points. 
Even with most of the freedom of the HOD model elim- 
inated, the modeling still reproduces the data with a re- 
duced statistic of ~ 0.4. We interpret this to indicate 
that the bootstrap errors estimated by lAllevato et al.l 
(|2011l ) are likely a conservative estimate of the true un- 
certainty. As a caveat to our modeling, however, we 
note that neglecting the covariance in the 2PCF mea- 
surements can also potentially lead to an overfitting of 
the data. 
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M [h-iMe] 

Fig. 3. — Comparison of the total mean occupation function 
of X-ray AGN as empirically measured with a group catalog by 
lAllevato etall ([2012; blue filled circles) to the 2PCF modeling in- 
ference from the current work (red solid line). The red shaded re- 
gion indicates the 68% confidence interval of our modeling (see the 
text). The black open circl es represent the empi rically-measured 
mean occupation function of lAUevato et af] 1)20121 ^ after correc ting 
for redshift and luminosity evolution of the X-ray AGN. See il4.2l 
for further discussion on the comparison between the two measure- 
ments. 



4.2. Comparison with Previous Work 

Our analysis suggests that the 68% confidence interval 
for the best-fit host halo mass of central X-ray AGN is 
12.90 < log[Mco„/(/i"^ Mq)] < 13.09. Using a simple 
HOD model consisting of only central AGN and charac- 
teriz ed by a delta functio n for the occupation distribu- 
tion, lAllevato et al.l (|2011| ) obtained a host halo mass of 
log[Mcon/(^"^ Mq)] = 12.9 7 ± 0.06, w hich is consistent 
with our result. Previously, iCoil et all ((2009) measured 
the CCF of X-ray AGN selected from the All- Wavelength 
Extended Groth Strip International Survey with a con- 
trol sample of galaxies and obtained a minimum host 
halo mass of x 10^^ Mq, which is also consistent 
with our result. The median redshift and luminosity of 
their AGN sample, 0.9 and 10*^^'^ erg s~^, r espectively, 
are similar those of the lAllevato et al.l (120111) sample. 

A r elated study was carried out bv iStarikova et ahl 
()2011|) . Using a sample of X-ray-selected AGN in the 
Chandra Bootes field (41 < log[Lx/(erg s"!)] < 45), 
they performed HOD modeling by comparing the cor- 
relation function projected parallel and perpendicular to 
the line of sight. For 0.17 < z < 3.0, their analysis pre- 
ferred a host halo mass scale > 4.1 x 10^^ h~^ Mq and 
found that AGN reside primarily in central galaxies. The 
derived host halo mass scale is consistent with our results 
but our satellite occupation is insufficiently constrained 
for a quant itative compari s on of satellite fractions. 

Recently, lAllevato et al.l ()2012D used a novel approach 
to directly measure the occupation function of X-ray 
AGN in groups and clusters. In Figure [3] we compare 



our total mean occupation function obtained from mod- 
eling the 2PCF (solid red line) to the direct measurement 
of lAllevato eTall (l2OT10(filled blue circles). The error- 
bars correspond to the la errors in mass and occupancy 
and the red shaded region represents the 68% confidence 
interval of our modeling. The best-fit values of our 
HOD parameters are as follows: \og[Mynin/{h~^ Mq)] — 
13.651^:^5, ^logM = 0.78l°:°9, log[Mi/(/i-i Mq)] = 

14.32l^ j^, and a = 2.59t? j^. It is evident from Figure[3] 
that the two complimentary measurements are consistent 
with each other, thus validating the HOD parameteriza- 
tion that we have used in the 2PCF modeling. As previ- 
ously mentioned, the model was developed using numer- 
ical simulations of lower luminosity AGN and has been 
applied to the current sample under the assumption of 
u niversality in th e gene ral AGN HOD properties. 

lAllevato et al.l (|2012| 1 attempt to correct their mea- 
surement of the mean occupation function for luminosity 
and redshift evolution. The measurements, represented 
by the black open circles in Figure [3l appear to disagree 
with our results at the 1.7a level. Physically, our {N{M)) 
can be interpreted as the mass-dependent AGN occupa- 
tion fraction at z ~ 1.2, under the assumption of weak 
luminosi ty and redshift evolution (see ij4.ip . We note 
that the lAllevato et al.l (2012) sample consists of AGN 
at z < 1.0 and wi th slightly low e r X-r ay luminosities 
than the sample of lAllevato et al.l (|201lD , whose 2PCF 
measurements we model in the current work. While it is 
possible that this disagreement in the occupation fraction 
could arise due to luminosity and / or redshift evolution of 
the X-ray AGN population, we emphasize that it could 
also arise simply due to an over- or under-correction for 
evolution. Conclusively identifying the source of this dis- 
crepancy will require studies of redshift and luminosity 
evolution with larger AGN samples. 

4.3. Implications for AGN Evolution 

From the measured clustering of AGN at different 
wavelengths and at different redshifts one can draw an 
evolutionary picture of AGN growth and connect it to 
the cosmic growth of structure (popularly known as 
'AGN co-evolution' in the literature). We now interpret 
our HOD result within this co-evolution p aradigm. We 
adopt the evolutionary picture presented in lHickox et al.l 
(|2009f l. which has been successful in interpreting sev- 
eral theoretical aspects of the co-evolution scenario, to 
assess our findings. W e will focus on the results of 
iRichardson et al.l (|2012D for the quasar HOD at z ~ 1.4 
and the results of the current paper for the X-ray AGN 
HOD at z ~ 1.2. Critically, the X-ray AGN we con- 
sider are subs tantially less bolom etrically luminous than 
quasars. The iHickox et al] (|2009[ ) picture has been eval- 
uated primarily on the basis of moderate-redshift obser- 
vations (z 0.5). We now extend this picture to a higher 
redshift {z ~ 1). 

In this picture, massive galaxies follow an evolution- 
ary sequence from a gas-rich star-forming quasar phase, 
to a spheroid with lower star formation and AGN ac- 

^ There is a typo graphical error i n the coefficients of the Lx ~ 
scaling relation in 'Alle yato et al.l I I2012I) . The correct coefHcionts 
((Leauthaud et al. 201(J) have been adopted in determining the 
masses of the X-ray groups and clusters in the paper (V. AUevato, 
private comm.). 
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Fig. 4. — Left panel: the median halo mass of quasars and X-ray AGN as a function of bolometric luminosity at 2: ~ 1, as labeled. The 
quasar result is obtained from Richardson ct al. (2012) and the X-ray AGN constraint is from the current work. Both values are obtained 
from HOD modeling of the 2PCF. This shows that there is a 2.5a difference between the host halo masses of quasars and X-ray AGN, with 
the latter residing in higher mass halos. This result is similar to t he evolution picture presented in |Hickox ct al. (2009) for AGN at z ~ 0.5, 
but now per formed at significantly h igher redshift (see also i|4.3l l. Right panel: the predicted total mean occupation function for quasars 
(dashed line: IRichardson et al.l [20121) compared to the predicted total mean occupation function of X-ray AGN (solid line; current work). 



tivity, to a red-and-dead galaxy wit h discrete radio out- 
bursts from the cen t ral AGN (e.g.. iCroton et al.l 120061 : 
iHopkins et all 120081: iHickox et all |2009 l ). The gas-rich 
quas ar phase could be possibly dri ven by major mergers 
(e-g.-lKauffmann fc Haehneltil2000l: [Hopkins ct al. 20063, 
|2008| ) or by secular instabilities (e.g.. iBower et al.i,2006f ). 
It is believed to occur at a typical halo mass scale of 
a few times 10^^ M© at almost all redshifts and may 
result from a natural evolution of the halo-galaxy mass 
relation, given that galaxy f ormation is most efficient at 
those halo mass scales (e.g.. iConrov fc Whiteli201 2') and 
powerful starbursts galaxies a re found in halos of similar 
mass (see IHickox et al.l 120121 and references therein) . In 
this phase the central black hole undergoes Eddington 
accretion and shines as a luminous quasar. 

From the quasar phase galaxies transition into the 
spheroid phase characterized by lower star formation and 
declining AGN activity. These spheroids reside in ha- 
los of mass ~ 10^'^ M©. Their central black holes un- 
dergo moderate accretion (0.001 — 0.1 Eddington) and 
shine as X-ray-bright AGN. The cause of this decline in 
AGN activity is not yet known, but a popular model 
involves a lack of cold gas due to shock heating in- 
side t he halo (e.g.. IBirnboim fc Dekell [20031 : IKeres et al. 
20091 ) or feedback froin the A GN itself (e.g.. lSilk fc ~ 
1998t lDi Matteo et al.ll2005[ ). At even later times, galax- 
ies enter the red-and-dead phase in which the central 
AGN becomes extremely radiatively inefficient. Con- 
tinued growth of the host halo completel y shock-heats 
the in tra-halo gas. From clustering studies IHickox et al.l 
(|2009f) demonstrated that X-ray AGN reside in dark mat- 
ter halos in the mass range of ^ 10^'^ Mq. Using 
a sample o f radio AGN, which have very low accretion 
efficiencies, IHickox et al.l (|2009f) also demonstrated that 



radio AGN are found in halos of mass ~ 3 x 10^"^ h ^ Mq 
or greater, c onsistent with the results of several other 
studies (e.g.. Wake et al.ll2008l: [Mandelbaum et al.ll2009[: 
iDonoso et al. 2010) 

IRichardson et all ()2012D found the median host halo 
mass of central quasars at z ~ 1.4 to be 4.11q4 x 
10^^ Mq, which is significantly {2.5a) lower than 
the corresponding host halo mass scale of X-ray-selected 
AGN near this redshift. The left panel of Figure [4| shows 
the repsective locations of quasars and X-ray AGN in 
(bolometric) luminosity-host halo mass space. It is ap- 
parent that, at ~ 1, X-ray AGN with low bolometric 
luminosities reside in halos of significantly higher mass 
than quasars with bolometric luminosities at least two or- 
ders of magnitude higher. Since, even at higher redshifts, 
higher luminosity AGN (i.e., quasars) reside in halos of 
lower mass compared to lower luminosity AGN, we con- 
firm that the correlation between AGN luminosity and 
host dark matter halo mass is relatively weak. Moreover, 
black holes indeed appear to undergo the quasar phase 
near the critical mass limit of ~ 10^^ Mq and generally 
become less active as the halo grow s in mass. Thus our 
detailed HOD modeling supports the lHickox et al.l ()2009f ) 
picture at high redshift. 

We note t hat t he "cartoon" picture presented in 
IHickox et al.l ()2009[ ) is valid only for central galaxies. The 
fact that our HOD modeling shows that AGN reside pri- 
marily in central galaxies validates interpreting the AGN 
population within this picture. An HOD analysis of the 
2PCF of high-redshift radio-selected AGN could further 
constrain this evolutionary picture. Some studies have 
suggested that, at high redshift, the CCF of radio-loud 
quasars with their radio-quiet counterparts does show 
more clustering strength than the auto-correlation func- 
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tion of radio quiet quasars (e.g.. lShen et al.|[2009D . which 
would indicate that radio-loud quasars reside in more 
massive halos. 

From the average duty cycle measurements of quasars 
and X-ray AGN we can infer a characteristic lifetime of 
these objects. The lifetime can be approximately esti- 
mated as to = fo X tu ([M artini fc Weinber^l2001l: [Crotonl 
l2009HHopkins fc Hernauist<,2009al ). where fa is the Hub- 
ble time and fo is the average duty cycle. For our pur- 
pose, the duty cycle can be estimated as the occupation 
fraction of AGN or quasars at the median redshift, which 
is shown in the right p anel of Figure |4l For quasars, 
(|Richardson et al.ll2012h obtained an average duty cycle 
of /o = 7.3t?j X 10-4 at z - 1.4. The Hubble time is 
roughly 6 Gyr at z ~ 1, yielding a characteristic quasar 
lifetime of ^ 4.4 Myr. Similarly, from the HOD model- 
ing results of the current work, we find a characteristic 
X-ray AGN lifetime of ~ 0.7 Gyr. Since these lifetimes 
invoke duty cycles calculated for the median host mass 
scale, we interpret these timescales as the mean lifetime 
of a "typical" quasar or AGN at z ^ 1. 

We can compare our estima tes of lifetime with 
the sem i-analytic predictions of IHopkins fc Hernguistl 
(|2009bD . Using an observed distrib ution of Eddington ra- 
tios and AGN model light curves, IHopkins fc Hernquist] 
()2009b[ ) found that a black hole spends different amounts 
of time at different stages of its evolutionary sequence. 
The lifetime at high Eddington ratios (> 0.1) is typi- 
cally 10 — 100 Myrs and at moderate Eddington ratios 
(0.001 — 0.1) is typically 0.5 — 1 Gyr. Our current sample 
of X-ray-selected A GN are likely to ha ve moderate Ed- 
dington ratios fe.g.. iHickox et al]|2009[ ) while the quasar 
population is representative of the high accretion stages 
of a black hole. Our estimated lifetimes, both for the 
highly accreting quasar phase and the moderately ac- 
creting X-ray phase, are in broad agreement with the 
IHopkins fc Hernguistl (|2009b( ) prediction. As we do not 
have a catalog of the masses of our black holes, this 
comparison is necessarily mostly qualitative in nature. 
It has been discussed in IHopkins fc Hernguisti (|2009bf ) 
that the lifetimes of black holes are a function of several 
parameters (e.g., black hole mass, luminosity, host halo 
mass). Hence a more accurate compa rison of lifetimes 
with the IHo pkins fc Hernguistl ()2009b h picture is possi- 
ble with black hole mass measurements of our quasar 
and X-ray AGN samples. Qualitatively, our results in- 
dependently support the idea that black holes spend sig- 
nificantly different times in different phases of their evo- 
lutionary sequence. On average, we find the lifetime of 
black holes in the moderate Eddington phase to be ^ 100 
times longer than that of the quasar phase. We sum- 
marize the key findings of our analysis in the following 
section. 

5. SUMMARY OF RESULTS 

We perform an HOD analysis of the projected 2PCF 
of X-ray-selected AGN at z ~ 1.2. We use a physically- 
motivated HOD model based on low-luminosity AGN in 
cosmological simulations, but which has also been found 
to provi de an excellent descript ion of the clustering of 
quasars ([Richardson et al.ll2012D and the mean occupa- 
tion empirically measured for X-ray AGN (l Allevato et al.l 
[2012i) . Our modeling yields a median mass for halos 



hosting central AGN of 1.02i[5J;^ x lO^^ /^-i ^nd 
an upper limit on the AGN satellite fraction of ^ 0.10, 
which are consistent with previous estimates. Our 2PCF 
data sample an insufficient range of the one-halo term to 
conclusively rule out a monotonically decreasing satel- 
lite occupation at high halo mass, but we find that our 
analysis strongly favors a p ositiv e slope. 

Since iRichardson et all ()2012D inferred a host halo 
mass of 4.1^0 4 x 10"'^^ Mq for central quasars at 
similar redshift, we show that X-ray AGN reside in halos 
of significantly {2.5a) higher mass compared to quasars 
at z '--^ 1. We also show that at z ~ 1, the average 
lifetime of the X-ray AGN phase is ^ 100 times longer 
than that of the quasar phase. The derived mass scales 
and lifetime estimates support the iHickox et al.l ()2009f ) 
picture of AGN evolution (established at low redshift), 
where black holes are believed to follow an evolutionary 
sequence from a high-Eddington quasar phase to a mod- 
erately accreting X-ray phase, to a radiatively inefficient 
radio phase. 

Based on our c urrent analysis , we pr oject the 2PCF 
measurements of lAllevato et al.l ()2011l : with appropri- 
ate rescaling of the error bars) to forecast the improve- 
ment in HOD modeling precision that can be achieved 
with future X-ray surveys. For a survey of equal area 
to XMM-COSMOS, we find that a sample of ~ 3, 000 
AGN is sufficient to constrain the HOD at the 20% level, 
while a sample size greater than 10,000 will enable char- 
acterization of the HOD to at least the 10% level (and 
possibly better, depending on survey systematics). How- 
ever, realizing this precision for the satellite occupation 
and its derived parameters will require measuring the 
2PCF to smaller pair separation scales. Based on the 
results of IRichardson et all (|2012D , we argue that this 
minimum scale is ~ 0.01 Mpc. This should be possi- 
ble with future X-ray surveys, whose larger volumes will 
a contain greater number of small-scale pairs. Our pro- 
jections show that ongoing and proposed X-ray surveys 
such as the Chandra COSMOS Legacy survey, a Chan- 
dra medium-area survey, and the eROSITA mission will 
easily obtain ^ 10%-level precision, provided that there 
is sufficient optical spectroscopic follow-up to obtain ac- 
curate redshifts. Accurately determining the AGN HOD 
will enable key physical parameters like the AGN satellite 
fraction and duty cycle to be inferred with unprecedented 
precision, making possible a definitive evaluation of the 
AGN co-evolution paradigm. 
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